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ABSTRACT 

Context. MWC 560 (V694 Mon) is one of the most enigmatic symbiotic system with a very active accretion-powered hot component. Such 
activity can be supported only by a luminous asymptotic giant branch star, i.e. a Mira or SR variable, with a high mass-loss rate. It is also a 
very unusual jet source because the jet axis lies practically parallel to the line of sight. 

Aims. The aims of our study are the determination of the evolutionary status of the cool component of MWC 560. 
Methods. Our methods involve analysis of near-IR JHKL and optical light curves. 

Results. The cool component of MWC 560 pulsates with a period of ~340 days, and it is probably a red SR variable on the thermally pulsing 
AGB. The high mass-loss rate expected for such a star is sufficient to power the observed activity of the hot companion. 



Key words, stars: binaries: symbiotic - stars: long-period variables - stars: mass loss ~ stars: individual: MWC 560 = V694 Mon - infrared: 
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1. Introduction 



Discovered by iMerrill & BurweUl (11943b as a star with strong 
emission lines, MWC 560 (V694 Mon) is an enigmatic sym- 
biotic system with a very active hot component. Its optical 
spectrum is always characterized by highly variable absorp- 
tion features, blue-shifted by 1000-6000 kms"'. These origi- 
nate from H I, He i, Ca ii and Fe ii, and are detached from the 
stationary, naiTow emission lines. The blue-shifted absorption 
features can be expl ained by a jet outflow along the line of sight 
dTomov et al.lll990l) . In 1990, MWC 560 underwent a 2 mag 
photometric outburst, accompanied by spectacular changes in 
these blue shifted absorption lines (e.g. Tomov et al. 199d 



Schmid et al ■1200 ll and references therein). The hot component 
is also a permanent source of rapid flickering with amplitudes 
in the range 0.2-0.7 m ag on time scales of 10-100 min (e.g. 
Gromadzki et al lllooel and references therein). 

Both the flickering and the 1990 outburst characteris- 
tics suggest that the hot component is predominantly pow- 
ered by unstable ac c retion, and its luminosity -100-1000 Lq 



dTomov et al.1 Il996t iMikolaiewski et al 



tively high accretion rates of M ^ 
200lh . 



19981) requires rela- 
10-^Moyr-i dSchmid et al.l 
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Doroshenko et al. I dl993h found a 19 30-day period in the 
mph historical light curve ( Luthardl]|l991 ) combined with more 
recent B-band photometry. They suggested an orbital origin of 
this period: if the orbit is eccentric, near periastron the accretion 
rate increases, and the hot component gets brighter. 

We know rather little about the cool component of 
MWC 560. However, the high activity of its companion 
requires a high mass-loss rate, consistent with an evolved 
AGB star rather than a normal red giant. Absorption fea- 
tures in the optical spectrum suggest an M3-M 4 spectral type 
dSanduleak & StephensonI ll973fclAUMil Il978l) . Infrared spec- 
tral classification mainly based on TiO ban ds, and the Call 
triple t indicates an M4-M5 giant dSzkodv e t al. 1990; BoppI 
19901; iThakar & Win3ll992l). Ba sed on the lack of VO bands 



at 1.05 jt/m Meier et al. classified the cool compo 

nent as an M5-M6 III giant. The most recent classification 
based on five TiO bands in the near infrared resulted in M5. 5 
and M6 (two different observations) (' Murset&Schmidll999h . 
[Frackowiak et al. (2003) found possible pulsations of red giant 
with f = 161 days, and concluded that it is an AGB star. 

This paper contains an analysis of near infrared (JHKL) 
and visual (AAVSO and ASAS) light curves with the objective 
of searching for periodic changes and establishing the nature of 
the cool giant. 
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Fig. 1. Light curves for MWC 560. From top to bottom: 
AAVSO+ASAS and JHKL light curves, and J - K colours, 
respectively. 



2. Observations 

JHKL broad-band photometry (1.25, 1.65, 2.2, 3.45 /jm) was 
obtained with the Mkll infrared pho tometer on th e 0.75 cm 
telescope at SAAO, Sutherland (see [Carter 1 19901 for details 
about the system). The measurements are good to +0.03 at 
JHK and +0.05 at L. The optical light curve consists of visual 
magnitude estimates collected by the American Association of 
Variable Stars Observers (AAVSO), and V-band photometry 
obtained in the fr ame of All Sky Automated Survey (ASAS) 
( PojmanskH 2002 ). The near-IR observations cover the period 
from November 1984 to February 2004 (92 points in / and 
K bands, 91 in // and 86 in L) whereas the optical photom- 
etry was collected between April 1990 and December 2001 - 
the AAVSO data (2051 points), and from November 2000 to 
February 2006 - the ASAS data (330 points), respectively. The 



near-IR data are listed in Table 1 and the light curves are pre- 
sented in Fig. 1. 

3. Period Analysis 

We analysed the near-IR, ASAS, and combined 
AAVSO+ASAS magnitudes by means of the program 
PERIOE rl ver 5.0, which use s the modified Lomb-Scargle 
d E 



method JPress & Rvbickil[l989l) . To remove long-term trends 
from the optical data, a second order polynomial was sub- 
tracted from the combined AAVSO+ASAS light curve, and a 
third order spline function was subtracted from the ASAS light 
curve. The resulting power spectra are shown in the right panel 
of Fig. 2, and results of our period analysis are summarized in 
Table 2. The periods given in the table were derived from the 
maxima of the peaks in the periodograms {f„l^), whereas their 
accuracy was estimated by calculating the half-size of a single 
frequency bin (A/), centred on the peak (fc) in a periodogram, 
and then converted to period units (AP = f^^ ■ Af). 

Power spectra of the composite AAVSO+ASAS light curve 
(Fig. 2, the uppermost, right panel) show a very pronounced 
peak at 1931'' (P ad) very close to the period reported by 
Doroshenko et alJ ([l993)- The second and third harmonics, 
~1931/2 and ~1931/3, respectively, are also present. We at- 
tribute this periodicity to repeated episodes of enhanced activ- 
ity of the hot component following the ephemeris: 

JD(max) = 2448 080+ 1931 x£ 

The light curve folded with this period is also show in Fig. 2 
(the uppermost, left panel). 

Power spectra corresponding to the JHKL light curves are 
very similar to each other (Fig. 2). In all spectra the strongest 
peak corresponds to a 339-day period (Ppui) which we take to 
be a real periodicity produced by pulsations of the cool compo- 
nent. There are also peaks around ~4500-5500 days (which is 
comparable to the time interval spanned by the data) connected 
with the long-term trend. Power spectra of the HK light curves 
show also peaks at ~ 1900'' days and ~ 310''. The former is 
very close to fact found in the optical light curve whereas the 
latter seems to be an annual alias of 1900-day peak rather than 
a real periodicity. In the power spectrum of the J light curve the 
~ 1 900-day and ~3 1 0-day periods are not seen because they are 
overwhelmed by the strong long term trend. These periods are 
absent in the power spectrum of the L light curve. The follow- 
ing ephemeris gives the phases of maxima in the JHKL bands: 

JD(max) = 2 445 960 + 339 x £ 

and the light curves folded with this ephemeris are shown in 
Fig. 2. (left panel). 

The power spectra of the visual photometry also show 
peaks at ~ 310''. These power spectra have poorer resolution 
that those for the near-IR data, so the peaks, especially that 
for the ASAS data, are relatively broad, and they may in fact 
be due to the 339'' pulsation period combined with one year 
alias of the 1931'' period. Such an interpretation is supported 



source of program is available on ^http://www.starlink. rl.ac. uk/\ 
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Table 1. IR photometry of MWC 560 



Table 1. continued. 



JD 


J 


H K 


L 


-2440000 








(day) 




(mag) 





6024.55 

6026.45 

6034.53 

6043.57 

6082.46 

6108.39 

6113.33 

6136.30 

6184.24 

6338.57 

6380.61 

6440.46 

6454.48 

6465.35 

7958.39 

7961.37 

7964.33 

7989.26 

8023.22 

8028.19 

8197.59 

8227.56 

8252.56 

8280.49 

8294.38 

8320.40 

8326.41 

8390.20 

8617.53 

8670.41 

8705.33 

8731.29 

8768.20 

8899.63 

8932.56 

8960.52 

8968.44 

8969.52 

8991.41 

9024.42 

9292.54 

9379.36 

9503.19 

9643.62 

9667.53 

9728.50 

10035.54 

10053.60 

10086.47 

10110.42 

10126.45 

10181.29 

10221.22 

10437.54 

10464.46 

10478.37 

10498.39 



6.50 
6.51 
6.48 
6.49 
6.54 
6.54 
6.56 
6.54 
6.49 
6.50 
6.49 
6.59 
6.56 
6.54 
6.43 
6.42 
6.43 
6.44 
6.42 
6.39 
6.47 
6.48 
6.48 
6.43 
6.43 
6.44 
6.47 
6.44 
6.44 
6.48 
6.49 
6.42 
6.45 
6.49 
6.47 
6.44 
6.42 
6.42 
6.44 
6.43 
6.47 
6.38 
6.52 
6.53 
6.51 
6.47 
6.48 
6.46 
6.50 
6.52 
6.52 
6.46 
6.45 
6.45 
6.49 
6.49 
6.48 



5.49 


5.12 




5.50 


5.14 


4.74 


5.49 


5.12 


4.79 


5.48 


5.13 


4.79 


5.52 


5.16 


4.84 


5.53 


5.18 


4.84 


5.56 


5.19 


4.86 


5.52 


5.16 


4.85 


5.44 


5.10 


4.90 



5.45 
5.54 
5.53 
5.50 
5.45 
5.44 
5.43 
5.42 
5.41 
5.42 
5.49 
5.49 
5.49 
5.44 
5.45 
5.44 
5.45 
5.43 
5.42 
5.41 
5.41 
5.37 
5.44 
5.48 
5.45 
5.43 
5.41 
5.43 
5.40 
5.41 
5.43 
5.36 
5.50 
5.52 
5.49 
5.44 
5.47 
5.44 
5.48 
5.49 
5.50 
5.46 
5.43 
5.44 
5.45 
5.46 
5.47 



5.12 
5.20 
5.17 
5.16 
5.11 
5.11 
5.10 
5.09 
5.06 
5.08 



5.08 
5.08 
5.07 
5.05 
5.03 
5.09 
5.14 
5.15 
5.10 
5.06 
5.09 
5.06 
5.06 
5.09 
5.01 
5.17 
5.15 
5.15 
5.10 
5.13 
5.10 
5.13 
5.14 
5.15 
5.11 
5.09 
5.10 
5.11 
5.12 
5.14 



4.73 
4.90 
4.89 

4.76 
4.80 
4.75 
4.78 

4.71 
4.84 
4.79 
4.82 
4.71 
4.74 
4.80 
4.79 
4.74 
4.77 
4.69 
4.70 
4.69 
4.74 
4.82 
4.85 
4.73 
4.71 
4.79 
4.70 
4.77 
4.79 
4.67 
4.80 
4.83 
4.82 
4.72 
4.81 
4.73 
4.79 
4.77 
4.87 
4.74 
4.72 
4.78 
4.79 
4.81 
4.86 



JD 


J 


H 


K 


L 


-2440000 










(day) 




(mag) 




10567.23 


6.50 




J. 14 


4.77 


10801.51 


6.54 


J Ay 


J. 14 


4.81 


10884.38 


6.47 




^ no 




11146.59 


6.48 


jAI 


^ 1 1 
J. 1 1 


4.86 


11187.48 


6.51 


jAq 


J.Li 


4.72 


11234.30 


6.46 


jAi 


J.Uo 




11479.54 






J. 14 




11505.51 


6.53 


Z AQ 

jAq 


J. 1 J 


4.77 


11511.48 


6.52 


^ Ai^ 

J AO 


J. Iz 


4.81 


11572.48 


6.52 


J. 40 


^ 1 1 
J. 1 1 




11599.35 


6.50 


^ Ai^ 


3. 1/ 


4.82 


11611.33 


6.52 


Z AQ 


3. 1/ 


4.78 


11617.34 


6.52 


^ Ad 

J Ay 


J. 1 J 


4.76 


11675.24 


6.50 


J. 40 


^ 1 1 
J. 1 1 


4.76 


11857.58 


6.54 


J.jl 


J. 1 / 


4.86 


11864.58 


6.55 


J.j4 


^ 1 O 

J. ly 


4.90 


11869.52 


6.56 




J. 10 


4.87 


11881.54 


6.54 


^ ^ 1 

J.jl 


J. to 


4.87 


11890.57 


6.55 




J. to 


4.85 


11927.47 


6.54 




J. to 


4.88 


11976.31 








4.85 


12239.58 


6.49 


J A/ 


J.Li 


4.75 


12255.57 


6.53 


^ AQ 

J Ay 


J. 1 J 


4.80 


12263.47 


6.51 


j.jyj 


J. Id 


4.75 


12284.50 


6.52 


z /in 
J Ay 


J. Li 


4.76 


12289.44 


6.51 


J. 4b 


J. Id 


4.81 


12319.40 


6.48 


jAj 


J. 12 


4.74 


12324.44 


6.49 


5.47 


5.10 


4.80 


12349.32 


6.48 


5.44 


5.10 


4.80 


12384.26 


6.49 


5.43 


5.09 


4.85 


12572.58 


6.59 


5.53 


5.17 


4.87 


12594.60 


6.56 


5.52 


5.15 


4.81 


12664.43 


6.52 


5.50 


5.15 


4.87 


12690.38 


6.52 


5.45 


5.09 


4.79 


12979.55 


6.49 


5.45 


5.12 


4.79 


13031.48 


6.48 


5.41 


5.06 


4.73 


13063.34 


6.41 


5.37 


5.03 


4.66 



by the presence in the ASAS power spectrum of another peak 
at ~ 166'*. This peak can be more accurately determined and 
it is due to the second harmonic of the pulsation period, which 
in this case would be 332*^. The ASAS light curve folded with 
this period is also shown in Fig. 2. 

In the power specti'um of the AAVSO-hASAS data there 
is also a strong peak at 747'' connected with a quasi-periodic 
oscillation which is predominantly visible during JD 2 450 500- 
2452 300 (1997-2002). However, we note that it is very close 
to 2 years, and it may be an artifact. 

4. Discussion and conclusions 

Our near-IR photometry began in 1984, with only a dozen mea- 
surements over a 1.5-year interval. The observations restarted 
in 1990, when MWC 560 underwent an outburst, reaching the 
brightest leve l in its whole photometric history, and ejected jets 
dTomov et al.lll99Q) . In the near-IR, it brightened by ~0. 1 mag 
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Table 2. Peaks in power spectra 



Frequency 


Period 


Power 


Remarks 


days"' 


days 


sigma unit 





J 








2.9479 10"^ 


339±4 


20.8 




H 








5.3276 10-" 


1877±126 


7.0 


Pact 


2.9479 10"' 


339±4 


19.8 




3.2321 10"' 


309±3 


10.6 


alias 


K 








5.3276 lO"'' 


1877±126 


7.2 


^act 


2.9479 10-3 


339±4 


20.2 




3.2676 10^* 


306±3 


12.2 


alias 


L 








2.9487 10"' 


339±4 


11.6 


Ppu[ 


A C A C 








3.2708 10"' 


306±12 


35.8 




6.0183 10"' 


166+4 


18.5 


f pui/2 


AAVSO+ASAS 








5.1789 10"-* 


1931±162 


285.6 


/'act 


9.4946 lO-* 


1053±48 


57.9 




1.3379 10"' 


747±24 


102.4 




3.1937 10-' 


313±4 


53.1 


Pas AS, alias 



in all, JHKL, bands. Then during the following 12 years, it 
slowly faded to the pre-outburst magnitudes observed in 1984. 
In 2002-2004, the star brightened again. The point scatter of the 
near-IR light curves was always of ~0. 1 mag in JHK bands, 
and ~0.2 mag in L band, respectively. In general, the near-IR 
light curves reflect the trend shown by the visual/V light curve, 
with additional maxima around JD 2 450 000 and JD 2 452 000. 
The J - K colour becomes redder as the brightness declines 
(Fig. 1) due to a decreasing contribution from the hot compo- 
nent. 

Our period analysis of the JHKL light curves revealed a 
339-day periodicity which can be attributed to radial pulsa- 
tions of the M giant. Although this period is close to that of 
Mira itself, and of many other Galactic Mira variables, the am- 
plitude of the pulsation, AK ~ 0.1 mag, is much lower than 
AK ^ 0.4 mag observed in Miras, including typical symbiotic 
Miras (IWhitelockl Il987h . We therefore classify it as an SRa 
variable (using the definition in the GCVS4) rather than as a 
Mira. This pulsation is hardly detectable in the visual light be- 
cause radiation in this range is dominated by the very active 
hot component. In fact, the scatter of points from the folded 
light curve in Fig. 2, AV ~ 0. 3 H- 0.6 mag is co mparable to 
the amplitude of the flickering (iGromadzki et al.ll2006i) which 
complicates the period analysis. The red giant pulsation period 
found in this stud y is more than twice the 161 -day pulsation 
period reported by iFrackowiak et al.l (l2003h . and it is the only 
periodicity present in all bands. 

The basic properties and evolutionary status of the semi- 
regular variables (SRVs) of type SRa and SRb were discussed 
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Fig. 2. Power spectra (right) and corresponding light curves 
folded with the strongest period (left). From the top to bot- 
tom: AAVSOh-ASAS fight curve folded witfi period Pact=1931 
days, ASAS ligfit curve folded with period Pasas=332 days, 
and JHKL ligfit curves all folded witfi period Ppui=339 days, 
respectively. 



in detail by iKerscfibaum & HronI d 19921 1 1994 Il996l) . In par- 
ticular, tfiey found tfiat tfie SRas appear as intermediate ob- 
jects between Miras and SRbs in all aspects, including peri- 
ods, amplitudes and mass-loss rates. Tfiey also concluded tfiat 
tfie SRas do not form a distinct class of variables, but are a 
mixture of 'intrinsic' Miras and SRbs. Tfie SRbs split into a 
'blue' group witfi P < 150 days and no indication of circum- 
stellar sfiells and a 'red' group witfi temperatures and mass- 
loss rates comparable to Miras, but periods about fialf as long. 
Tfiey suggested tfi at tfie 'red' and 'Mira' SRb s are tfiermally 
pulsing AGB-stars (IKerscfibaum & Hronlll992l) . Tfie persistent 
and relatively long period places tfie cool component of MWC 
560 among tfie 'Mira' SRVs. Tfiese differ from normal Miras 
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only in their smaller pulsation amplitudes. The IRAS [12//m]- 
[25 fim]-0.64 locates MWC 560 in the period-IRAS colour di- 
agram of Kerschbaum & Hron ( 1 992) in the region occupied by 
86 % of the Miras, whereas most SRbs are outside this region. 
The average, reddening corrected (/sb-v = 0.15; Schmid et al. 
2001) neai--IR colours {H - K}o = 0.30, {J - H}q = 1.00, 
K - [12 yum]=1.25, and the spectral type M5.5-6 are also con- 
sistent with an O-rich SRV, although the colours are unlike 
those of a Mira. Finally, if the 1931 day periodicity is orbital 
then the two stars of this symbiotic system are close to each 
other (see below), and the presence of a nearby companion 
may influence the pulsation characteristics, possibly reducing 
the pulsation amplitude. We suggest that the evidence supports 
the view that the red component of MWC 560 is on the TP 
AGB -phase. 

The SRV component of MWC 560 could be the source of a 
strong stellar wind and thus support the observed high activity 
of the hot component. If we estimate the mass of the red giant 
at IMq and that of the white dwarf companion at O.5M0 then, 
assuming that the 5.3 year period is orbital, the separation of 
the two stars will be 3.55 AU and the radius of the Roche lobe 



for the M giant in a circular orbit will be 1 .5 AU. ISchmid et al 



(1200 ih argued that an accretion rate up to a few x 10~^ M0 yr 
is necessary to account for the observed hot component lumi- 
nosity which is comparable to the win d efliciency found for 



typical Miras and SRVs. For example, lOlofsson et al.l (l2002h 
determined average mass loss rate of 2 x 10"^ Mq yr"', with a 
maximum value of 8 x 10"^ Mq yr"' for a sample of M-type 
irregular and semi-regular variables. So, in the case of MWC 
560, a significant fraction of the mass lost in the cool giant 
wind must be accreted by the hot companion. This is possible 
only if the binary components of MWC 560 interact via 'wind 
Roch e lobe' overflow occur s instead of a spherically symmetric 
wind (|Podsiadlowskill2007l) . 

One of the most intriguing features of this hot component 
activity is its highly periodic character The 1931 -day period- 
icity found in the visual AAVSO-i-ASAS light curve is essen- 
tially identical with t he 1930-dav periodicity detected in the 
OTpg/B light curve by iDoroshenko et al.l (Il993h . We note that 
this periodicity remained in-phase over a century, and the most 
natural explanation would be an orbital motion. The only prob- 
lem with the orbital interpretation is that the orbital inclination 
of MWC 560 seems to be extremely low (as indicated by the jet 
axis practically aligned with the line of sight, and the lack of or- 
bitally related radial velocity changes) which excludes any ge- 
ometrical effects, such as eclipses, illumination, etc. Although 
one can argue that the jet axis can be inclined with respect 
to the binary orbital plane, the relatively large amplitude of 
this variability, Ampg ~ 1 and AV ~ 0.5 mag, respectively, 
points to another mechanism(s). One possibility is an eccen- 
tric orbit and enhanced accretion rate near periastron passage 



causin g a brightening of the hot component ( IDoroshenko et al 



19931) . Then we can even speculate that if the eccentricity is 
high enough, ~ 0.5 or more, a Roche lobe overflow may occur 
near the periastron. 

Another possibility is enhanced mass loss due to periodic 
changes in the SRV environment. We note that the 1931 day 
period equals roughly six times the pulsation period. A peri- 



odically enhanced mass loss might be caused by some kind 
of interplay between the cool component pulsation and lay- 
ered dust formation shells. Such effects have been shown by 
time-dependent hydrodynamic simulation of C-rich Mira envi- 
ronments (Winters et al. 1999; Winters & Le Bertre 2001). In 
particular, these simulations have been able to reproduce vari- 
ations of the mass-loss rates on different time scales. 
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